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B and 5 C mesons with NRQCD and Clover actions 
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We present preliminary results from our study of the heavy-light spectrum and decay constants. For the heavy 
quark, we use NRQCD at various masses around and above the b quark mass. For the first time, the heavy 
quark action and the heavy-light current consistently include corrections at second order in the non-relativistic 
expansion, as well as the leading finite a corrections. The light quarks are simulated using a tadpole-improved 
Clover action at various masses in the strange and c quark region. 



1. INTRODUCTION 

Heavy-light hadrons are of topical interest in 
theoretical and experimental particle physics. 
The B meson decay constant fs is required, to- 
gether with the bag parameter Bb , for a determi- 
nation of the poorly known elements of the CKM 
matrix, in particular, of the phase leading to CP- 
violation in nature. Comparison with experimen- 
tal results on the heavy-light spectrum available 
currently and in the near future provides an im- 
portant test of the reliablity of the approxima- 
tions used in these calculations. 

Using NRQCD, it is possible to simulate the b 
quark directly, without having to use extrapola- 
tions in the heavy quark mass. Earlier results 
for /VM, obtained using NRQCD consistent 
through 0(l/rag) where itiq is the bare heavy 

quark mass (see e.g. [0j2|), showed that fy/M has 
a large slope as a function of 1/M. Such a behav- 
ior would be an indication that higher order terms 
are important. We address this issue here. 



sate for the leading 0(a) error in the discretiza- 
tion of the temporal derivative. Furthermore, we 
include the leading relativistic correction, which 
is expected to be the largest 0(l/(mg) 3 ) correc- 
tion. This action is also consistent through 0(v 4 ) 
in the heavy quark velocity in heavy-heavy sys- 
tems. 

We calculate the heavy quark propagator G t 
on a timeslice t using the explicitly time-reversal 
invariant evolution equation: 
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The free heavy quark Hamiltonian Hq contains 
only the kinetic energy term 



Ho 



AW 



(2) 



and the other terms are similar to 



2. THIS SIMULATION 

In this study, the heavy quark action and 
heavy-light currents are consistently improved 
through 0(1/ (nig) 2 ) at tree level. We also in- 
clude the leading 0(a 2 ) lattice spacing correction 
to the lattice laplacian, and a term to compen- 
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Here A, A^l and A' 4 ' denote the tadpole- 
improved symmetric lattice discretization of the 
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spatial covariant derivative, the laplacian and the 
lattice-artifact fourth derivative described in ||]. 
The temporal gauge link, the magnetic field B 
and the electric field E are also tadpole-improved 
(i.e. each link is divided by uq = 0.878). E + and 
E~ denote the traceless forward and backward 
2-leaf clover electric fields, respectively. 

The heavy light currents (J) are also improved 
to the same order: 

J = Jo + ^ + 8^ J2 - (4) 
where, at tree level, 

J = qTQ, J 1 =qT(i-A)Q, 
J 2 =qr(A^> + a- B-2i^-E^jQ, (5) 

q being the light spinor and Q the nonrelativistic 
heavy quark spinor whose upper two components 
are zero. For axial vector currents, T = 7570, and 
for vector currents, r = 7^. Here E denotes the 
standard 4-leaf clover electric field. 

In addition, we have calculated the matrix ele- 
ments of the other operators with the right lattice 
symmetry and up to the mass dimension of J 2 - 
These are expected to mix under renormalization 
with the tree level operators. As the renormaliza- 
tion constants (and hence the mixing terms) are 
currently unknown, we do not discuss them any 
further. We simulate heavy quarks with masses 
am° Q = 1.6,2.0,2.7 with n = 2; 4.0,7.0,10.0 with 
n = 1, and in the static approximation. The three 
lower mass values are in the region of the b quark, 
and the others are used to study the extrapolation 
of the decay constant to the static limit. The pre- 
liminary results presented here are based on an 
ensemble of 57 quenched Wilson configurations, 
fixed to Landau and then to Coulomb gauge, on 
a 16 3 x 48 lattice at (3 = 6/g 2 = 6.0. To increase 
statistics we also include the time reversed ver- 
sions of these configurations. For light and charm 
quarks we use a tadpole improved clover action. 
We use two K values (0.119 and 0.126) around 
the charm region and three (0.1369, 0.1375 and 
0.13808) around the strange quark mass. 

The heavy quark smearing functions are 
hydrogen-like wave functions with zero or one 
node, or delta functions. All combinations of 



Table 1 



Pion masses with degenerate clover quarks. 



K 


0.119 


0.126 


0.1369 


0.1375 


0.13808 


am p 


1.526(2) 
1.552(2) 


1.162(2) 
1.200(2) 


0.424(3) 
0.554(6) 


0.363(3) 
0.51(1) 


0.297(4) 
0.47(1) 



those smearing functions at source and sink are 
used. The light quark propagators are always 
smeared at the source with a Gaussian smearing 
function with zero or one node, whereas at the 
sink we use either a smearing identical to that 
at the source, or a delta function. The results 
presented here are obtained with the zero-node 
smearing of the light propagator at the source. 

3. THE SPECTRUM 

3.1. Light spectroscopy 

In Table @, we list the masses of pseudoscalar 
and vector mesons with degenerate clover quarks. 
These masses are determined from the fall-off of 
the zero-momentum meson correlator and thus 
for the heaviest quarks, where one expects sig- 
nificant discretization effects in the dispersion re- 
lation, they should only be taken as rough es- 
timates. Using the method described in [0, we 
find the critical k to be k c — 0.13924(4) (in 
agreement with the UKQCD estimate ||), the 
value of the isospin averaged light quark kappa 
K[ = 0.13917(4), and the strange k to lie between 
0.1372(4) (determined from K* or $ mesons), and 
0.13755(19) (from K mesons). The inverse lattice 
spacing from m p is a -1 = 1.9(1) GeV. 

3.2. Quarkonia 

For each of the heavy quark masses used for 
the heavy-light simulation, we also calculate the 
S and P wave states of the corresponding heavy- 
heavy mesons, which is computationally quite in- 
expensive. Results will be presented elsewhere. 

3.3. Heavy-Light: B and B c 

For the B and B c spectrum, we study 1 Sq, 3 Si, 
3 P , 3 Pi, 3 P 2 and 1 P 1 states. The data on 1 P 1 - 
3 Pi mixing will be reported elsewhere. For the 
3 p2 states, operators belonging to both the E and 
the T 2 representations of the cubic group are used. 
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3.3.1. Meson masses 

To calculate the mass of the heavy light 1 S'o 
mesons, we first determine the energy difference 
E s im(p 2 ) — -Esim(O) by fitting the ratio of non-zero 
and zero momentum meson correlators to an ex- 
ponential. This is then related to the mass M2 by 

E sim ( P 2 ) - E sim (0) = yW-Ma - M 2 • (6) 

The values for M2, determined using p 2 = 1 and 
p 2 = 2 (in units of (27r/16a) 2 ) with k = 0.1369, 
are presented in the second and third columns 
of Table ||. For comparison, the fourth column 
lists results obtained assuming that the non- 
relativistic dispersion relation holds for p 2 = 1 , 

E aim (p 2 ) - E sila (0) = ■ (7) 

Alternatively, the mass shift A = M2 — E snn 
can be calculated perturbatively. We use A cal- 
culated at lowest order in a s and 1/toq ||, but 
for the action used in ||. E s - lia (0) is the rate of 
exponential fall-off of the zero momentum corre- 
lator. These M2 are given in the fifth column of 
Table |. 

Additionally, one can try to determine A from 
the degenerate heavy-heavy S state mesons using 

A=±A HH = ±(M 2 -E sim ) HH . (8) 

Determining M2 with the nonrelativistic dis- 
persion relation (ffl) for the heavy-heavy mesons, 
and E^~ from the fall-off of the heavy meson 
correlator, we obtain the last column of Table ||. 

We find that all three methods agree for the 
lighter masses within errors. For the heaviest 
masses, the deviations are ~ 10 — 20%. These 
are larger than the expected uncertainty of ~ 0.1 
in A pert from higher orders in perturbation the- 
ory, and indicate the presence of other systematic 
errors. Thus, in the remainder of this paper, we 
use M2 from the first method. 

3.3.2. Hyperfine splitting 

We extract the hyperfine splitting AE directly 
from the exponential fall-off of the ratio of the 
3 S*i and the 1 S'o correlators. For the three light 
k values we don't observe any variation of AE 



Table 2 



1 5'o meson masses at n — 0.1369. 



am Q 


p 2 = l 


p 2 = 2 


aM 2 
NR 


pert 
m Q 




1.6 


2.16(2) 


2.13(3) 


2.18(2) 


2.13 


2.20(3) 


2.0 


2.55(2) 


2.52(3) 


2.56(3) 


2.52 


2.60(5) 


2.7 


3.23(4) 


3.19(5) 


3.24(4) 


3.22 


3.33(7) 


4.0 


4.49(7) 


4.40(10) 


4.51(8) 


4.60 


4.7(1) 


7.0 


7.75(11) 


7.27(15) 


7.79(8) 


7.21 


8.4(2) 


10.0 


11.0(2) 


10.2(3) 


11.02(16) 


10.06 


13.1(2) 



Table 3 

Hyperfine splitting in lattice units. 



am Q 


0(l/{m%) 2 ) 


0(l/m£) 




k = 0.119 


k = 0.1369 


k = 0.137 


1.6 


0.0164(3) 


0.0198(10) 




1.71 






0.020(1) 


2.0 


0.0141(3) 


0.0165(9) 


0.017(1) 


2.7 


0.0116(2) 


0.0127(7) 




4.0 


0.00865(15) 


0.0086(6) 


0.0091(8) 


7.0 


0.00548(12) 


0.0046(4) 




8.0 






0.0042(7) 


10.0 


0.00865(15) 


0.0030(4) 





with k, so we present our data without any chi- 
ral extrapolation. Even in the charm region, the 
dependence of AE on n is small. In Table | we 
list the results at K = 0.119 and k = 0.1369, as 
well as the results obtained in yj] with an action 
correct to 0(l/rag). For the B mesons, we ex- 
pect the results for the two actions to agree, since 
the higher order terms in Eq. || do not affect the 
hyperfine splitting. The data, shown in Fig. [I], 
substantiate this. Thus the predictions for the 
hyperfine splittings in jl| remain unchanged. 

3.3.3. P states 

A typical effective mass plot for a P wave state 
with a heavy quark around the b mass and a light 
quark around the charm mass is shown in Fig. ||. 
For lighter k values, the signal is only slightly 
noisier. The data presented in Table ^ and shown 
in Fig. ||, indicate that the separations between 
the P states are qualitatively as expected. Note, 
however, that as the mixing between 1 P\ and 3 Pi 
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Figure 1. Hyperfine splitting as a function of the 
mass of the 1 Sq meson presented in column two 
of Table ^, with a linear fit to the data. Circles 
are new results, whereas the squares are from § 
using an action correct to only 0(1/toq). 



Table 4 

P - 1 S splittings for m° Q = 1.6, k = 0.119. 



0.20- 



0.15- 



state 


Lattice Rep. 


AS 


a Po 


Ai 


0.178(13) 


'Pi 


Ti 


0.203(12) 


3 Pi 




0.193(8) 


3 P 2 


£ 


0.204(9) 


3 P 2 


P 2 


0.200(12) 



i p 



3 Pi 3 P 2 (P) 3 P 2 (T) 



Figure 3. P — 1 Sq splittings lattice units (n 
0.119, rriQ — 1.6). Errors are purely statistical. 
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Figure 2. 

k = 0.119, = 1.6. Both the light and the 
heavy quarks use no-node smearing functions at 
the source, whereas at the sink only the light 
quark is smeared. 



has been ignored, thus the 1 Pi — 1 S'o and 3 Pi — 1 Sq 
splittings are not individually reliable. Also, even 
though the statistical errors are of the order of the 
3 P2 — 3 Po splitting, a study of the ratio of those 
correlators indicates that both lattice representa- 
tions of the 3 P 2 are equal and heavier than the 
3 Po. 

4. DECAY CONSTANTS 

4.1. Analysis 

The decay constants are related to the one 
meson production amplitudes of the local pseu- 
doscalar or vector current Jl- 



fVM = V2A L = yf2/M{0\J L \B) 



(9) 



To determine Al, we fit to the smeared-smeared 
and smeared-local 1 S'o and 3 Si correlators with 
different operator and smearing combinations at 
the source and sink to the exponential form: 



(OsnkO S r C ) = A Qsrc: O snk exp(-P 



src,snk 



t). (10) 
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The measured matrix of the amplitudes is then 
factorized as Ao arc ,o snk — Ao STC Ao snk by min- 
imizing the overall % 2 . The matrix element for 
the improved current, specified in Eq. |[ is then 
constructed by simply adding the various terms. 
To get the physical decay constants we extrap- 
olate linearly in 1/(2k) to k = k; and k = k s . 
The entire procedure is carried out in a single- 
elimination jackknife loop to estimate the errors 
in the final quantities. 

The signal for the effective amplitudes of the 
correlators show clear plateaus except for the 
static case. In Fig. ||, we show an example of 
the effective amplitude of the largest 0(l/(mg) 2 ) 
current correction to the decay constant. 




-100 1 1 1 1 1 1 1 1 1 1 — L 

10 20 

t 



Figure 4. Effective amplitude for the correlator 
of a 1 Sq current correction at itiq = 1.6 and k = 
0.1369. Both the quarks use smearing functions 
with no nodes at the source, and are unsmeared 
at the sink. 



4.2. Results 

We present a study of the effect of improving 
the action and the operators to 0(l/(mg) 2 ). It 
is important to recall that all coefficients have 



Table 5 

The contributions to the decay constants from 
0(1/toq) current correction (second column) and 
the three 0(1/(toq) 2 ) current corrections for k = 
0.1369. All numbers in the table have been mul- 
tiplied by 100 for convenience. 





(j-D 


a ■ B 


a-E 


D 2 


m° Q 


2m% 




4(m° Q ) 2 


«K) 2 


k = 0.1369 


1.6 
2.0 
2.7 
4.0 
7.0 
10.0 


-5.40(19) 
-4.59(17) 
-3.59(14) 
-2.62(12) 
-1.97(07) 
-1.58(06) 


0.82(04) 
0.49(03) 
0.25(02) 
0.10(01) 
0.035(3) 
0.018(1) 


1.16(05) 
0.75(03) 
0.41(02) 
0.19(01) 
0.060(5) 
0.031(3) 


-1.23(04) 
-0.86(03) 
-0.53(02) 
-0.27(01) 
-0.134(5) 
-0.080(3) 


k = 0.119 


1.6 

2.0 
2.7 
4.0 
7.0 
10.0 


-8.10(17) 
-7.18(16) 
-6.00(16) 
-4.67(14) 
-3.04(08) 
-2.28(07) 


1.28(03) 
0.78(02) 
0.41(01) 
0.18(01) 
0.055(2) 
0.025(1) 


1.60(03) 
1.07(02) 
0.62(02) 
0.30(01) 
0.106(3) 
0.054(2) 


-3.27(07) 
-2.44(05) 
-1.61(04) 
-0.90(02) 
-0.39(01) 
-0.231(7) 



tadpole-improved tree-level values with no per- 
turbative corrections. Thus, the following results 
should be considered qualitative and preliminary. 

The data in Table [| show that in the mass 
region of the B meson, (jtoq ~ 2.0), the indi- 
vidual current corrections arising at 0(l/(mg) 2 ) 
are ~ 20% of the 0(l/mg) corrections. How- 
ever for the pseudoscalar decay constant the three 
0(1/ (toq) 2 ) terms add up to a result close to zero 
because they have different signs. (For the vec- 
tor decay constant these three terms come with 
the same sign). The data for the different heavy 
quarks is shown in Fig. ||. 

The effects of the 0(l/(wig) 2 ) corrections in 
the hamiltonian are shown in Fig. ^| where we 
compare the chirally extrapolated results against 
the 0(l/mg) estimates given in 0. We find that 
the 0(l/(mq,) 2 ) data indicates a smaller slope, 
though the difference is marginally significant. 
We are in the process of calculating the pertur- 
bative corrections to the coefficients in order to 
extract the physical value of the decay matrix el- 
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Figure 5. Comparison of the chirally extrapo- 
lated values of fy/M with various current cor- 
rections: circles include all corrections through 
0(l/(mg) 2 ), crosses include up to 0(l/mg) cor- 
rections, whereas squares are the leading order re- 
sults. The fits shown are uncorrelated quadratic 
fits: broken lines show the error on the fit to the 
circles. 
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1 / aM, 



Table 6 

(fy/M)ps at various levels of the current correc- 
tions (no corrections, correct at 0(1/ 'ttiq) and at 
0(l/(mg) 2 )). The last column shows results of a 
calculation where both the Hamiltonian and cur- 
rent are correct to only 0(1/toq). 



m Q 


0(1) 


o(K)- 1 ) 


0((m Q )-) 


111 


1.6 


0.178(09) 


0.153(9) 


0.157(9) 




1.71 








0.150(11) 


2.0 


0.180(10) 


0.158(9) 


0.160(9) 


0.154(12) 


2.7 


0.178(11) 


0.161(11) 


0.161(11) 




4.0 


0.181(12) 


0.169(12) 


0.168(12) 


0.182(13) 


7.0 


0.189(18) 


0.178(18) 


0.178(18) 




8.0 








0.21(2) 


10.0 


0.216(25) 


0.208(24) 


0.208(25) 




oo 


0.249(9) 


0.249(9) 


0.249(9) 


0.25(4) 



ements and their slope as a function of 1/M. 
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Figure 6. The chirally extrapolated fvM cal- 
culated consistently at 0(l/mg) from |jj (data: 
squares, fit: dotted line) and at 0(l/(rag) 2 ) 
(data: circles, fit: solid line). The error in the 
0(l/(mg) 2 ) points is shown by a dashed line. 



